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Abstract. Cell polarization plays a central role in the development of complex 
organisms. It has been recently shown that cell polarization may follow from 
the proximity to a phase separation instability in a bistable network of chemical 
reactions. An example which has been thoroughly studied is the formation of 
signaling domains during eukaryotic chemotaxis. In this case, the process of 
domain growth may be described by the use of a constrained time-dependent 
Landau-Ginzburg equation, admitting scale-invariant solutions d la Lifshitz and 
Slyozov. The constraint results here from a mechanism of fast cycling of molecules 
between a cytosolic, inactive state and a membrane-bound, active state, which 
dynamically tunes the chemical potential for membrane binding to a value 
corresponding to the coexistence of different phases on the cell membrane. We 
provide here a universal description of this process both in the presence and 
absence of a gradient in the external activation field. Universal power laws are 
derived for the time needed for the cell to polarize in a chemotactic gradient, and 
for the value of the smallest detectable gradient. We also describe a concrete 
realization of our scheme based on the analysis of available biochemical and 
biophysical data. 



PACS numbers: 64.60.My, 64.60.Qb, 87.16.Xa, 87.17.Jj, 82.39.Rt, 82.40.Np 

1. Introduction 

Biophysical processes of cell polarization have attracted large interest in recent times. 
It has been observed that intriguing similarities exist in the polarization of such diverse 
biological systems as cells of the immune system, social amoebas, budding yeast, and 
amphibian eggs [38| . This suggests that cell polarization may be a highly universal 
phenomenon. 

One of the best studied examples of the role of biochemical cell membrane 
polarization in eukaryotic cells is chemotaxis. Chemotaxis is the ability of cells to 
sense spatial gradients of attractant factors, governing the development of all superior 
organisms. Eukaryotic cells are endowed with an extremely sensible chemical compass 
allowing them to orient toward sources of soluble chemical signals. This mechanism 
is the result of billion years of evolution, and multicellular organisms would not exist 
without it. Slight gradients in the external signals produced by the environment 
induce the formation of oriented domains of signaling molecules on the cell membrane 
surface. Afterwards, these signaling domains induce differentiated polymerization of 
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the cell cytoskeleton in their proximity, inducing the formation of a growing head and 
a retracting tail, and eventually directed motion towards the attractant source. 

It has been suggested in the biological literature that domains of signaling 
molecules are self-organized structures |25| . In this paper we confirm that this 
expectation may be substantiated by the use of statistical mechanical methods, leading 
to the prediction that universal features typical of coarsening processes in phase- 
ordering systems should be observable in polarizing cells. We also describe here a 
concrete realization of our scheme in the process of eukaryotic chemotaxis, based on 
the analysis of available biochemical and biophysical data. Part of the results presented 
here have been briefly reported in a previous letter [12] . 

2. Cell polarity 

Stochastic reaction-diffusion systems are a natural paradigm for describing in physical 
terms the biochemical processes taking place in the living cell, since the cytosol and cell 
membrane arc inherently diffusive environments^. Although active transport processes 
also take place in the cell, they regard mainly vesicles, organelles and large multiprotein 
complexes, while smaller cell constituents move diffusively. Thermal agitation and the 
intrinsic stochasticity in the advancement of chemical reactions provide natural sources 
of noise. 

Most reactions in the cellular environment would be very slow if they were not 
favored by the action of catalysts. Small numbers of enzymatic molecules (10^-10^ 
per cell) control the speed of chemical reactions involving much larger numbers of 
substrate molecules (10^-10^ per cell.) Often, the substrate concentration in its turn 
controls the catalyst activity, so that the response of the system becomes nonlinear. 
Most biochemically relevant reactions involve enzyme-substrate couples and arc part 
of networks of interconnected autocatalytic reactions. 

Nonlincaritics allow in principle the system to realize several stable biochemical 
phases, characterized by different concentrations of chemical factors [36]. Transitions 
between different phases in reaction-diffusion systems have been observed in purely 
physical settings, such as the adsorption and reaction of gases on catalytic surfaces 
[27l [40] . Recently, it has been shown that a similar process of nonequilibrium phase 
separation may be at the heart of directional sensing in higher eukaryotes [TTl [T^] . 

In eukaryotic directional sensing cells exposed to shallow gradients of 
external attractant factors polarize accumulating the phospholipidic signaling 
molecule phosphatidylinositol trisphosphate (PIPS) and the PIP3-producing enzyme 
phosphatidylinositol 3-kinase (PI3K) on the cell membrane side exposed to the highest 
attractant concentrations, while phosphatidylinositol bisphosphate (PIP2) and the 
PIP2-producing enzyme phosphatase and tensin homolog (PTEN) accumulate on the 
complementary side [24] (see [Appendix A| for a more abstract description of the 
relative roles of these signaling molecules.) 

Accurate quantitative experiments [331 128j performed by exposing Dictyostelium 
cells to controlled attractant gradients showed that uniform concentrations of external 
attractant factor induce a predominant, uniform concentration of PIP3 and PI3K on 
the cell membrane, and do not immediately result in cell polarization and motion. 
However, slight gradients in the distribution of the attractant factor induce the 
formation of two complementary domains, one rich in PIP3 and PI3K, and one rich 

X For general facts regarding cell biology wo refer to Rof. [2]. 
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in PIP2 and PTEN, in times of the order of a few minutes. This early breaking of the 
spherical symmetry of the cell membrane induces cell polarization and motion |24j . 
Uniformly stimulated cells observed over longer timescales (of the order of 1 hour) are 
seen to polarize stochastically and move in random directions. 

Numerical simulations of a stochastic reaction-diffusion model of the process 
suggest that both the early, large amplification of slight attractant gradients and 
the separate phenomenon of late, random polarization under uniform stimulation are 
explained by the proximity of the system to a spontaneous phase separation driven 
by non- linear autocatalytic interactions [llj . In this framework, cell polarization is 
the final result of a nucleation process by which domains rich in PIP3 and PI3K are 
created in a sea rich in PIP2 and PTEN, or vice versa, depending on initial conditions 
and activation patterns. The polarization process is accomplished when pure PIP2 
and PIP3 rich domains grow to sizes comparable to the size of the cell. Gradient 
activation patterns strongly influence the kinetic of domain growth and coalescence, 
taking advantage of the underlying phase-separation instability. This way, the peculiar 
reaction-diffusion dynamics taking place on the surface of the cell membrane works as 
a powerful amplifier of slight anisotropics in the distribution of the external chemical 
signal. 

In this statistical mechanical point of view, random and gradient-driven 
polarization appear as two faces of the same coin, in good agreement with some of the 
existing biological intuition [38] . 

To better understand the process of spontaneous and gradient-driven cell 
polarization from a physical point of view it is convenient to describe the corresponding 
signaling network in abstract terms, i.e. forgetting about the particular nature of the 
molecules involved and considering only the general structure of the network. This 
approach has the potential to provide a unified description of polarization phenomena 
in distant biological systems. 

In our abstract signaling network (Fig. [1]) a system or receptors transduces an 
external distribution of chemical attractant into an internal distribution of activated 
enzymes h, which catalyze the switch of a signaling molecule between two states, that 
we denote here as ip~ and A counteracting enzyme u transforms the if'^ state 
back into ip^ . The molecule ip^ in turn activates w, thus realizing a positive feedback 
loop. 

The signaling molecules ip^ , ip^ are permanently bound to the cell surface S 
and perform diffusive motions on it, while the u enzymes are free to shuttle between 
the cytosolic reservoir and the membrane. In a more complete description we should 
consider that also the h enzymes are shuttling from the cytosol to the membrane 
[TT| [12] . Here however for simplicity we represent with h only the receptor-bound 
fraction, which we identify with the external activation field. The diffusivity of u 
enzymes in the cytosol is much higher than the diffusivity of (p^ , (p^ molecules on 
the cell membrane, therefore membrane-bound u enzymes may be assumed to be in 
approximate equilibrium with the ip'^,ip~ concentration field. This fact leaves only 
the surface molecule concentration as relevant dynamic variables. Moreover, 

since the ip'^,ip~ molecules may only be converted into each other, we are left with 
only one relevant degree of freedom, their difference tp = ip'^ — tp~ . 

The model of Fig. [T]was initially introduced to describe chemotactic polarization 
in higher eukaryotes [11] [T^. In that case, we identify (p~ and (p~^ with PIP2 and 
PIP3, u with activated PTEN, and h with activated PI3K. 

Recently, it has been proposed that polarization of budding yeast (a lower 
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Figure 1. Local structure of a prototypical signaling network for cell polarization. 
Here </3+ and tp^ represent the local concentration of distinct signaling molecules, 
or distinct states of the same molecule, which are converted into each other by 
the couple of counteracting enzymes h and u. The u enzymes are activated by 
, resulting in an amplification loop. The surface distribution of h enzymes is 
assumed to simply mirror an external distribution of soluble chemical attraetant. 
The signaling molecules , i/p^ are permanently bound to the cell membrane and 
perform diffusive motions on it, while the u enzymes are free to shuttle between 
the cytosolic reservoir and the membrane. 



eukaryote) may be the result of an amplifying feedback loop similar to the one 
described in Fig. [1] [4]. In our language ip~ and ip^ represent there the activated and 
unactivated states of the Cdc42 small GTPase (see Appendix Ap , while u would be 
identified with the activating factor Cdc24. The model of Ref. [1] lacks a counteracting 
enzyme playing the role of h in the scheme of Fig. [1] and is therefore not bistable. For 
this reason, it can reproduce only stochastic, intermittent polarization, as is observed 
at the border of the bistability region in the case of chemotactic polarization [Til [H] ■ 
However, in a recent work [34] a counteracting Cdc42 deactivating factor that could 
play the role of h has been described. This suggests that polarization of budding 
yeast cells may be driven by a bistable potential allowing the realization of stable 
polarization, similarly to the case of chemotactic polarization of higher eukaryotes. 



3. Macroscopic description of cell polarization 

Cell polarization is a macroscopic effect, emerging from the stochastic dynamics of 
a network of chemical reactions taking place in occasion of the random encounters 
of specific signaling molecules which perform diffusive motions and shuttle between 
the cell cytosol and membrane [26l [181 EH] ■ A large amount of information has been 
collected in recent years about the biochemical aspects of cell polarization in higher 
[261 [m [381 [241 [231 [351121] and lower [3ll [Ml [2Q] eukaryotes. However, available data 
cannot be considered yet complete or quantitative to a satisfactory degree. This kind of 
situation is typical of present efforts to derive macroscopic aspects of cell behavior from 
noisy and yet poorly quantitative data about the relevant microscopic interactions. 
It is therefore extremely important that a sensible macroscopic description of cell 
polarization can be given, starting only from the knowledge of a few robust properties 
of the biophysical system. 

In this Section we develop such a description. In the next Section wc show how 
known examples of cell polarization fit in our general scheme. 

Let us start here by assuming that we have knowledge only about the following 
robust properties of the cell polarization process: 
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1. Single component order parameter: The state of the system may be effectively 
described in terms of the configurations of a single-component concentration field 
(f describing the distribution of a set of signaling molecules on the cell membrane S. 

2. Bistability: The underlying chemical reaction networks allows the realization of 
distinct, locally stable chemical phases. 

3. Self-tuning: A global feedback mechanism controls the metastability degree ip of 
the system and drives it towards a state of phase coexistence. 

4. Non-conserved field: There are no local constraints on the values assumed by the 
field (p. 

The present set of properties stems from the abstraction of known properties of 
eukaryotic polarization (sec also the next Section). In particular. Property 4. is 
the consequence of the fast diffusion of u enzymes across the cytosolic volume |12j . 
(It is worth mentioning here that our framework would still hold, although with a few 
differences, also in the case that Property 4. be substituted by a local conservation 
condition.) 

Property 1. implies that the evolution of the state of the system can be described 
by a single stochastic reaction-diffusion equation. Studies of non-equilibrium statistical 
mechanics have shown that a few classes of nonlinear stochastic equations may emerge 
from the coarse-graining of microscopic dissipative dynamical systems, depending on 
general properties, such as the number of field components and the presence, or 
absence, of local conservation laws |15[[7]. 

Property ^. leads us to select the time- dependent Landau- Ginzburg model 

dMr.t)^-'-p^ + Eir,t) (1) 
Sip{r,t) 

(or model A in the classification of Rcf. [15 ) where 



dr (2) 



is an effective free energy functional, h is an external activation field, 13 is a diffusion 
constant, V^^h is an effective potential, and S is a noise term taking into account the 
effect of thermal agitation and chemical reaction noise. 

Property 2. implies that the effective potential V^^h has two potential wells, 
corresponding to a couple of distinct, stable chemical phases ip+ and </3-[§l. The kinetic 
advantage of transforming a region of phase into a region of if - phase is measured 
by the metastability degree 

The polarized state corresponds to the stable coexistence of the </5+ and ip- phases in 
complementary regions of the cell membrane. 

Property 3. implies that -0 is an integral functional of the field configuration, 
going to zero for large times under stationary conditions. A reasonable analiticity 
assumption then leads to the following system of equations, describing the dynamics 
of cell polarization in the presence of a stationary external activation field: 

dMr, t) = Z3VV(r, t) - ^ [^(r, t)] + S(r, t) (3) 

^(t) (X / tp (r, t) dr — / (y9(r, oo)dr, t ^ oo (4) 
Js Js 

§ We are using a slightly different notation to distinguish the values assumed by the ip field 

from the names of the concentration fields ip~^ , ip~ of signaling molecules. 
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It is possible to derive a concrete realization of the scheme described in the previous 
Section in the case of the signaling network of Fig. [1] by using the law of mass 
action, the quasistationary approximation for enzymatic kinetics, and the limit of fast 
cytosolic diffusion (see [Appendix B\ . In this case, the state of the system can be 
described by the single concentration field (p = (p^ — (p~ , thus giving Property 1 of the 
previous Section. The field is not constrained by a local conservation law because 
ip'^ molecules can be freely converted into Lp~ molecules and back on any point of the 
cell surface. This fact corresponds to Property 4- 

The evolution of the p field is described by the equation 



dtp = DS/'^tp - kce^tKassf- 



2K + C + P 2K + C-P 

where / = Uftee is the volume concentration of free cytosolic u enzymes (which is 
approximately uniform as a consequence of fast cytosolic diffusion), ft, is a surface 
activation field, ifass is the association constant of u enzymes to (p~ signaling 
molecules, /ccat is a catalytic rate, K is a saturation (Michaelis-Menten) constant. 

The corresponding effective potential has the form Vf^h{p) ~ fViif) + hV2{p) 
(see [Appendix B[ ) . The metastability degree ■(/; is therefore a function of h and /. If 
h = h(r, t) is not uniform {e.g. if the cell is exposed to a chemical activation gradient) 
tp takes on different values in different points of the membrane surface. Wc consider 
however for the moment the simplest case where the activation field is uniform in 
space and constant in time. 

A simple analysis shows ( [Appendix Bp that there are regions of parameter values 
such that y^./i is bistable, with two potential wells ip>^ and p- corresponding to stable 
phases respectively rich in the p^ and p^ signaling molecules. Thus Property 2 is 
verified. 
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Figure 2. Typical form of the effective potential V^^fi{ip) (red curves) in different 
parameter regions (c/. supplementary text of Ref. Region I: one equilibrium 

point ip- < c; Region II: one equilibrium point ip+ = c; Region Ilia: one stable 
equilibrium 73„ < c and one metastable equilibrium = c. Region IIIj,: same 
as Ilia, but now ip^ = c is stable and ip— is metastable. The phase coexistence 
curve separating Ilia from III;, is defined by the condition V^^hiv—) = Vilj,h{v+)- 
Arrows show the direction of the dynamic drift towards the phase coexistence 
curve. 
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In the present problem, the volume concentration / of free enzymes varies in 
time (but not in space) . More information about its values can be gotten in the limit 
(realized for small membrane diffusivitics and large times) when the interface between 
the and the phase is much smaller than the typical domain size, allowing us to 
use the so-called thin wall approximation. Then, the value of / is simply linked (see 
[Appendix B[ ) to the area covered by the (p- phase: 



showing that the metastability degree is proportional to the excess fraction of free 
cytosolic u enzymes with respect to their value at equilibrium. The presence of this 
global feedback mechanism corresponds to Property 3. 

The present situation is closely reminiscent of the decay of the uniform, metastable 
state of a supersaturated solution with the formation of precipitate grains |19| . In that 
case, the metastability degree is proportional to the excess solute concentration with 
respect to its equilibrium value. The main difference with the present case is that in 
the case of precipitation, the density field is locally constrained by the law of particle 
conservation [7] , and its evolution is described by Model B of Ref. [15] , instead than 
by model A. 

5. Phase separation kinetics 

In polarization experiments cells are exposed to uniform or gradient distributions 
of attractant factors and polarize either spontaneously, or in the direction of the 
attractant gradients [38j . The properties of the model free energy described in the 
previous Section and numerical simulations of a model system suggest that the 
introduction of an external attractant distribution moves the system in a region of 
bistability, where the uniform phase realized at initial time becomes metastable and 
germs of a new phase are nucleated. Depending on the way the system is prepared at 
initial time, the metastable phase can be either a ip^ rich or a p~ rich phase. 

The process of decay of a metastable state in physical systems described by 
systems of equations similar to ([31 S]) has been extensively studied in the framework 
of the theory of first-order phase transitions [HI [7]- The process passes through 
successive stages of nucleation, coarsening, and coalescence (Fig. [3]). In the first stage, 
approximately circular germs of the new, stabler phase are produced in the sea of the 
metastable phase by random fluctuations, or by the presence of nucleation centers. 
In the second stage, a process of coarsening is observed, where larger domains of the 
new phase grow at the expense of smaller ones, the average size of domains grow, and 
the average number of domains decreases. In a finite system, the process is concluded 
when a state of phase coexistence is reached. In this final state, the two phases are in 
equilibrium and are polarized in two large complementary domains. 

For our purposes, a detailed knowledge of the initial, nucleation stagefj] is not 
necessary, as long as its characteristic time t^ is so fast that a large number of germs 
of the new phase is nucleated all over the cell surface, well before the coarsening stage 
start sH. 

II And therefore of the precise characteristics of the noise term H which is its driving force. 
^ The converse case, where to is the largest timcscale of the problem and polarization is the result 
of the rare nucleation of a solitary domain, cannot provide a mechanism of gradient sensing which is 
at the same time insensitive to the uniform component of the attractant field, and highly sensitive 
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To understand the subsequent, coarsening state we have to focus on the laws by 
which the domains of the new phase either grow or shrink. 

We consider here the case when the new phase is a minority phase, so that we can 
restrict our consideration to approximately circular domains, which are dominating 
because they minimize the linear tension between the two phases. For simplicity, we 
shall also restrict to domains which are small enough that membrane curvature may 
be neglected. 

An approximate equation for the growth of a circular domain of size r may 
be derived from ^ in the thin wall approximation. Inserting the approximate 
propagating solution (^(R, = (j){R — r{t)) (Fig. [4^) for the radial domain profile 
in ([3]) and integrating over S we get 

^•^v>a[0] ^ _ ^ _^ ^/ (5) 
dr dr 

where 

Q = y / (0')'dR 
^ s 

is a dissipation function |17j and ^' is a noise term. 

For a circular domain of radius r, Q ~ ■firrr'^, where 

7= (0')'di?= J ^2V^.M{^)/Dd(j) 

is a kinetic coefficient [3. On the other hand the effective free energy for a circular 
domain or radius r is [7]: 

J-,p^h = Zttot — nr'^'ip (6) 
where a = is a linear tension. 

From (0 in]) we get the following approximate equation for the growth of a circular 
domain of size r: 

-yr = h ^ (7) 

r 

where ^ is a noise term. 

Eq. ([7]) shows that domains smaller than the critical radius 
a 

rr = — 

are mainly dissolved by diffusion, while germs with r > r^ mainly survive and grow 
because of the overall gain in free energy (Fig. Wp)- 

to its gradient component. Indeed, the nucleation of a single domain could provide a mechanism of 
gradient sensing only if the gradient would induce significantly different domain nucleation rates in 
different points of the cell membrane. But in that case, also variations in the uniform component 
of the attractant field would produce large variations in the typical polarization times, while the 
converse has been reported. 
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Figure 4. a: Radial profile of a growing (r > rc) and a shrinking (r < Vc) 
circular domain of the phase in the sea of the phase. The and the 
ip- phase are separated by a diffusive interface of thickness <5. b: Qualitative 
graph of the effective free energy of a circular domain of size r. Domains larger 
than r > Tc tend to grow because the energetic gain due to the surface term 
nr^i(> overcompensate the loss due to a longer interface. The converse happens 
for domains with r < rc. 



During the nucleation stage the noise term produces a population of germs of the 
new phase of size close to 

rg ^ Tc (5 

in a characteristic time to. For domains with r > ?'c the noise term in ([7]) may be 
neglected and domain growth is an almost deterministic process. 

It is interesting to estimate ro ^ 6 in terms of observable parameters. The 
thickness 5 can be estimated as 

where b is the potential barrier separating the two phases [7]. The height of the 
potential barrier may in its turn be estimated dimensionally from (jB.14p as 6 kcnthc, 
giving 

Using realistic parameter values (D ^ 1//to^/s, fccat ^ c/h ^ 10) we get 



6. The coarsening stage 



When domains of the new phase occupy an appreciable fraction of the membrane 
surface S a coarsening stage sets on. Domain growth makes the degree of metastability 
tjj decrease and renders further growth of the new phase more and more difficult. The 
critical radius rc grows with time, so that domains that earlier had size larger than rc 
become undercritical and shrink, and larger domains grow at the expense of smaller 
ones. In a large system rc soon becomes the main length scale in the problem, leading 
to the appearance of a scaling distribution of domains of size r. 

The population of coarsening domains of size r can be described in terms of the 
size distribution function n{r, t), such that n{r, t)Ar is the average number of domains 
with size comprised between r and r + Ar, and the total number of domains at time 
t is given by 

/•oo 

N{t)= / n{r,t)dr 
Jo 
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The time evolution of n{r, t) implied by ([7]) is dcseribed by a standard Fokkcr-Planck 
equation |36| . If we restrict our consideration to supercritical domains we can neglect 
the diffusive part of the Fokker-Planck equation since for them the noise term ^ is 
negligible. This means that the stochastic nature of the problem enters mainly in the 
formation of the initial distribution of germ sizes n{r^to), while for r > rc the time 
evolution of n{r, t) is dictated by the deterministic part of Thus, we are left with 
the following kinetic equation: 



dn{r, t) d 
^ dt ^d^ 



[m-^)nir,t)] =0 (9) 



Eq. © contains the unknown function ip{t), and is therefore not closed. We 
obtain a closed system by complementing ^ with the asymptotic law 



oo 

2 



^{t)(xAoo- 7rr'^n(r,t)dr (10) 



obtained from ^ in the thin wall approximation. Here 

/•oo 

Aoo = / 7rr^n(r, oo)dr 
Jo 



is the area occupied by the new phase at equilibrium. 

For large times a scaling distribution of domain sizes can be found explicitly 
( [Appendix"^ and Fig. [5]): 

n{r,t)dr = p{r/rc)d{r/rc), ^p{t) = — 

rc = rc{t) ^ ro{tltoY''^ (11) 

where 

where ro is the characteristic domain size at the beginning of the coarsening stage and 
C~0.11. 

The total number of domains decreases in time due to the evaporation of small 
domains. Using the explicit solution pTl [T2|) . we easily find: 

C C Aqq I rfj 



Mp) = 77^^exp(-;^), to = ^, (12) 



N{t) = I n{r,t)dr 







r? t 



Similarly, it is possible to compute explicitly the value of the average domain size, 
which is found to coincide exactly with the critical radius: 

{r) = rc 

7. Spontaneous and gradient- induced polarization 



The coarsening theory exposed in the previous Section allows to deduce a simple 
scaling law for the time needed for spontaneous cell polarization. 
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12 3 4 



r/rc 

Figure 5. Time evolution of tlic sclfsimilar domain size distribution n(r, t) 
{t/to = 1,2,3,4). 

If the cell has size R, the growth of domains according to pTjl comes to a stop 
at the time when the average patch size (r) becomes of the order of the cell size R. 
From (|lip we get 

- to {R/rof 

At the end of the process the cell is polarized in a random direction. The actual 
direction of polarization is the result of the initial random unbalance in the germ 
distribution. 

The typical time for random polarization is of the order of 10'^ s [12]. Together 
with the estimate (|5]) this gives to ^ 10 s. 

Let us now consider the case where a source of external attractant is present at 
some distance from the cell, in such a way that a gradient of external attractant is 
created by diffusion close to the cell surface (Fig. [S]). 




source 



Figure 6. Geometry of stimulation of a cell by an external source of attractant. 

The inhomogencity in the distribution of attractant induces a similarly 
inhomogcneous distribution of activated enzymes h. This way, the degree of 
metastability tp takes on different values on different points of the cell surface. 

If the cell membrane has a nearly spherical form and a radius R much smaller than 
the characteristic scale of the attractant distribution, and if the gradient component of 
the activation field is small with respect to the background component on the scale R, 
the metastability degree ip at the beginning of the coarsening process may be written 
as the sum of a uniform component ip and a small space-dependent perturbation: 

ip + Sip with 6ip = —eipo cos 9 

where ipo is the value of the uniform component at the beginning of the coarsening 
process and e is the relative gradient on the scale R. The perturbation modifies the 
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equation of domain growth ([7]) as follows: 

jr = ip eipo cos + ^ 



(13) 



where 9 is an azimuthal angle defined in Fig. [51 

The uniform component tp varies in time together with the (approximately) 
uniform concentration of u molecules in the cell volume. On the other hand, the 
perturbation Sip is constant in time, but not uniform in space, being proportional to 
the external attractant distribution. 

As long as etpQ <C ipi the effect of the perturbation is negligible, so domain growth 
proceeds according to the law pT|) and the uniform component ip decays as 

In a large cell there is a crossover time t^ when the perturbation becomes of the 
same order of the uniform component: 

ip{t^) = etpo 

Using the scaling law (fTTj) we get 

to 

c2 



After domain growth enters in a new stage, where the growth becomes anisotropic. 
Domains in the front and back of the cell get different average sizes (Fig. [7]). 



nucleation 



coaxsemng 




coalescence 



t > u 



Figure 7. The crossover time separates an initial stage of isotropic coarsening 
from a final stage when domains evaporate from the back of the cell and condense 
in the front. 



Indeed, for t > t^ the leading term in ([T^ is the perturbation eipo cos 9, implying 
that in the region closer to the source of the perturbation (cos 9 > 0) the (p- phase 
evaporates, and in the region away from the source (eos6' < 0) it condenses. At the end 
of the process, complete polarization is realized (Fig. \7^. In this final stage domains 
grow approximately linearly in time, thus the total time t'^ to reach polarization is still 
a quantity of order (using definition p4p from the next Section it can be estimated 
as i {1 + e/eth)tt). 

The above scheme is valid as soon as the initial nucleation time to is significantly 
smaller than t^, an assumption which is compatible with the observation of real |25] 
and numerical |11| experiments. 



8. Gradient sensitivity 

The second stage of domain evolution described in the previous Section occurs only 
if > t^. Otherwise, the presence of a gradient of attractant becomes irrelevant and 
only the stage of isotropic domain growth actually occurs. This condition implies that 
a smallest detectable gradient exists, such that directional sensing is impossible below 
it. The threshold value eth for e is found by the eonditionfe = t*. Since the product 
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Figure 8. Dynamic mechanism of gradient sensing. In the presence of a gradient 
(red) domain sizes grow initially as t^^'^. After the crossover time te they grow 
linearly in time and the cell polarizes in the direction of the gradient (see also 
Fig. (TJi. In the absence of a gradient (blue) the t^^'^ growth goes on until the cell 
polarizes in a random direction. 



tjjrc is a time- independent constant, we can simply compare its value at initial and 
final time when e = eth, obtaining that the threshold detectable gradient is 

..h = ^ (14) 

Using the estimates from Sections [5] and [3 and the typical value R ^ 10 /im, we get 
e ~ 10%, a value which is compatible with the observations [55] , 

An interesting speculation is that the bound (|14[) may explain why spatial 
directional sensing was developed only in large eukaryotic cells and not in smaller 
prokaryotes, whose directional sensing mechanisms rely instead on the measurement 
of temporal variations in concentration gradients [3] . By solving ()14p in terms of the 
size R we get the following bound for the size of a cell which may be able to sense a 
relative gradient e: 

e 

Our bound goes in the same direction as the size criterion formulated in [5], but it's 
independent of it, since the criterion of Rcf. [5] is based on estimates of signal-to- noise 
ratios, while our bound stems from the intrinsic properties of polarization dynamics. 



9. External fluctuations 



One may wonder whether a cell may become polarized by transient gradients produced 
by a spontaneous fluctuation in the external distribution of attractant molecules, 
or fluctuations in receptor-ligand binding, as has been suggested in the literature 
[18] . Since eukaryotic cells typically carry 10^-10^ receptors for attractant factors, 
one expects spontaneous fluctuations in the fraction of activated receptors to be of 
the order of 10^, a value which is comparable to observed anisotropy thresholds. 
However, to actually produce directed polarization the fluctuation should sustain itself 
for several minutes, i.e. for a time comparable to the characteristic polarization time 
(such as ). Such an event has very low probability of being observed since the 
correlation time of the fluctuations determined by attractant diffusion at the cell 
scale and the characteristic times of receptor-ligand kinetics are much less than the 
polarization time. Indeed, the diffusion time is ~ 1 s at the typical cell size 10 /im, and 
the characteristic times of receptor-ligand kinetics are also ~ 1 s (see online supporting 
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information to Rcf. [33]). Therefore, the direction of ceU polarization in the case of a 
homogeneous distribution of attractant can only be determined by the inhomogeneity 
in the initial distribution of the positions of PIP2-rich germs produced by thermal 
fluctuations. 

10. Conclusions 

By using standard statistical mechanical methods we have shown that the dynamics 
of signaling domains in cell polarization is independent on the nature of the signaling 
molecules and the values of kinetic rate constants, as long as some very general 
conditions are met: 

a) Timescale separation allows to describe the polarization process in terms of a single 
concentration field of signaling molecules on the cell membrancQ- 

b) The underlying chemical reaction network is bistable. 

c) A global feedback mechanism drives the system towards phase coexistence. 

d) The cell is sufficiently larger than the size of nucleating germs of the new phase. 

These conditions allow the cell to work as a detector of slight gradients of external 
stimulation gradients. 

The property of universality arising from our analysis cannot be underestimated. 
Presently, several efforts are made to understand the dynamical behavior of 
living beings starting from microscopic informations provided by molecular biology. 
However, these informations are mostly incomplete and poorly quantitative, and 
theories that depend in a sensitive way on them are likely to be of little utility. But 
if some behavior happens to be universal, a consistent physical theory of it may be 
built, which can be compared to experiments. 

The universal properties of cell polarization emerge from properties of domain 
growth which have been extensively studied in first-order phase transitions [7]. The 
similarity of the two problems follows from the fact that fast degrees of freedom 
of chemical kinetics are in approximate equilibrium with slower degrees of freedom, 
which can be described by means of an effective free energy functional. It is worth 
observing that in the biological system studied here, there is no direct interaction 
between signaling molecules, similar to the one observed in solid state system such as 
binary alloys, but only an effective interaction mediated by enzyme activity, binding, 
unbinding and diffusion processes. 

Our theoretical scheme allows to shed light on some non trivial questions, such 
as the mechanism of directional sensing and the effect of random fluctuations of the 
medium on the polarization process. Random polarization appears as the result of 
the intrinsic stochasticity of the process of domain nucleation and not of random 
fluctuations of the medium. Random and gradient-induced polarization appear as the 
two sides of a same coin. Our scheme provides an explanation of why spatial directional 
sensing is not observed in the small prokaryotic cells, and provides asymptotic 
estimates for polarization times and threshold detectable gradients. 

An important component of our picture is the existence of a global coupling of 
the degree of mctastability to the state of the system [121 [IHl ■ The constrained phase- 
ordering dynamics tunes the system towards phase coexistence, similarly to what 

+ Wo should consider adding hero the condition that the concentration field is not locally constrained 
by a conservation law. However, also the converse case of a locally conserved field can be treated in 
a similar way without substantially changing the present scheme. 
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happens in the case of a precipitating supersaturated solution. The global control 
allowing self-tuning to phase-coexistence is realized by shuttling of enzymes from the 
cytosol to the cell membrane and backwards. 

Some of the features that we have observed in cell polarization have been 
considered in previous works, such as the fact that equations of the form ([T]) are 
relevant for the description of systems of bistable chemical reactions [29l [36] . and 
that global couplings in activator-inhibitor reaction-diffusion systems may lead to the 
formation of stable spatiotemporal patterns [HI [30] . The peculiar properties of this 
kind of systems have led to the use of the term of excitable or active media. Using 
this same language, we can say that the cell membrane acts as an active medium 
responding to the stimulation with the formation of domains of a new phase. Our 
work proposes that directional sensing results from the peculiar, universal features of 
the phase-ordering dynamics of these domains. 

From a biological point of view, the universality of the polarization process 
allows the cell to behave in a robust, predictable way, independent on microscopic 
peculiarities such as the precise values of reaction rates and diffusion constants. 

We first proposed that chemotactic cell polarization may result from the simple 
ingredients of bistability induced by a positive local feedback loop in a signaling 
network and global control induced by shuttling of enzymes between the cytosol 
and the membrane in our previous works [TTJ [TH [H] ■ Other authors have proposed 
similar models, either independently [32] or subsequently [21] (a review of models of 
chemotactic polarization can be found in Ref. [E]). Some of these models try to 
take into account computationally the interactions of a large numbers of chemical 
factors, while retaining the essential role of a feedback loop as generator of a phase- 
separation instability. However, most of the reaction rates that should be provided 
to perform such computations are known with very poor accuracy. Our framework 
suggests however that such a detailed description may be not necessary, as long as 
properties [a|,...(d| are met. 

Aspects of the bistable mechanism of eukaryotic polarization firstly introduced in 
Ref. [TT] (supporting material) have been considered in recent papers [51 [22] as relevant 
to polarization phenomena. A similar mechanism, out of the bistability region, has 
been proposed to explain intermittent polarization in budding yeast [4] . These works 
suggest that the combination of bistability and global control [Tl] [12] is providing a 
useful paradigm for the understanding of cell polarization phenomena. 

Acknowledgments We thank Guido Serini for many inspiring discussions. This 
research was supported in part by the National Science Foundation under Grant No. 
NSF PHY05-51164. 

Appendix A. Lattice gas description of cell polarization 

The signaling molecules PIP2 and PIP3 are different phosphorylation states of 
the phosphatidylinositol molecule, i.e., they carry a different number of phosphate 
groups attached (2 and 3, respectively). Enzymes which catalyze phosphorylation 
of their substrate, i.e. the addition of a phosphate group, are called kinases, while 
dephosphorylating enzymes are called phosphatases . 

It is natural to visualize the state of a chemical system such as the one described 
in Fig. [Tjin terms of two families of classical spins on a twodimensional lattice, taking 
on values -1 (PIP2, PTEN), (an empty site), +1 (PIP3, PI3K) [10]. Taking into 
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account fast cytosolic diffusion, the enzyme family becomes slaved to the substrate 
family [TU]. 

In this lattice-gas description the existence of a cytosolic enzymatic reservoir 
exchanging enzymes with the cell membrane is represented by a chemical potential 
for enzyme creation and destruction (actually, adsorption and desorption to/from the 
cell membrane), globally coupled to the lattice configuration [TO] , 

The PIP2 and PIPS molecules constitute approximately 1% of the total number 
of membrane phospholipids, and the number of PI3K and PTEN enzymes are at least 
one order of magnitude lower, thus, both the substrate and the enzyme population 
should be thought as diluted gases. 

Two-state (or multistate) molecules such as PIP2 and PIPS are all but an 
exception in cell biology. Another example is given by small GTPases, such as the 
Cdc42 molecule involved in the polarization of budding yeast, which can be found 
either in the activated GTP state or in the deactivated GDP state. The switch 
between the two phosphorylation states is catalyzed by a couple of activating (GEF) 
and deactivating (GAP) enzymes [2j. 



Appendix B. Mean-field equations for eukaryotic polarization 

We derive here mean-field equations for eukaryotic polarization using standard 
methods of chemical kinetics, including Michaelis-Menten saturation terms for the 
enzymatic componcnts0. We make use of the fact that the diffusivity i^voi of u 
enzymes in the cytosol is much faster than the diffusivity D oi ip molecules on the cell 
membrane: this fact allows to considerably reduce the number of dynamical degrees 
of freedom. 

We describe the macroscopic state of the cell using surface concentration fields of 
membrane-bound molecules (Fig. [1]) and the volume concentration field / = Wfrco of 
free u enzymes. 

The chemical kinetic equations for the signaling network of eukaryotic polarization 

are: 

dtip+ ^ Dy^if+ - fceat , , + fccatTT^ (B.l) 

K + (p+ K + (p- 

dtf' = i?VV- + fccat J"^^ , ~ fccaty^T dtu (B.2) 

dtU fcass/v" - kdissU (B.S) 

dtf =V-(i?™iV/), (B.4) 
They must be complemented with the boundary condition 

df 

J = D,,i^^dtu (B.S) 
an 

where d/dn is the derivative along the outward normal to the membrane surface S. 
Condition (|B.5|) expresses the fact that the fiux of u enzymes leaving the cytosolic 
volume equals the flux of enzymes being bound to the cell membrane. 

For simplicity, we consider here identical catalytic, association and dissociation 
rates (fccat, fcass, fcdiss) ^ud Michaelis-Menten constants K for the (p^ (p~ and 

— > ip'^ processes. This is compatible with existing information about these 

* Michaelis-Menten saturation terms arise from timescale separation in enzymatic kinetics, which 
allows to make use of a quasi-stationary approximation [5]. 
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processes, suggesting that reaction rates differ by factors of order 1 [TT] and allows 
to easily study the equations analytically. 

Typical values for surface and cytosolic diffusivity are D 1/im^/s, I?voi ~ 
lQ^m?/s [11| . Typical values for rate constants are: fccat ~ ^diss ^ I'S"^, fcass ~ 
0.05s~^nM~^; for the total number of (p^ and Lp~ molecules, and the total number of 
u and h enzymes: ~ 10^, Ny, Nh lO'^-lO^. Observe that Nu/N^ < 1. 

The usual definition of macroscopic fields such as u is as follows. For each point 
r in space we choose a volume v centered in r, containing n{v) molecules, and we 
compute concentrations as u(r) = liuiy^Q n(v) / v . This implies that the number of 
molecules of the relevant chemical factors is so large that v can be chosen much 
smaller than the size of the system, but large enough that the resulting field i^(r) is 
approximately continuous. This hypothesis is not always acceptable, since enzymatic 
molecules are present in the cell in very small numbers. We shall therefore assume 
that real concentrations are described as the sum of an average part u, described 
by mean field equations of the kind (jB.lHB.5)) . and a fluctuating part 5u taking into 
account both the discrete character of the concentration field and thermal disorder. 
The fluctuations du due to random adsorption and desorption processes are at the 
origin of the noise term S in ([3]) (see [Appendix C[ ). 

Since enzyme diffusion in the cytosol is faster than phospholipidic diffusion on 
the membrane, during the characteristic times of the dynamics of membrane-bound 
factors, /(r, t) relaxes to the approximately uniform value 

/W = ./o--^ / u^{r,t)dr (B.6) 



V „ ^ 

where /g = N^/V , while u relaxes to the local equilibrium value 

U = Kassfv' (B.7) 
where i^ass = fcass/fcdiss- 

On the other hand, by summing jRU and (|R2| we get 

dt (<^+ + (^-) = DV^{p+ + (^-) - dtu^ (B.8) 

Since Nu/N^p <C 1 we neglect the term dtu. Then. (|B.8p shows that the sum 
c = (p^ + if" tends to be approximately uniform and constant in time. 

By subtracting (jB.l[) and (|B.2p and introducing the difference concentration field 
if = ip~^ — ip~ we get 

dtp = DV p- fccat :r^r- ■ + fccat TTTT- B.9) 

2K + c + p 2A + c — p 

and using the local equilibrium condition (|B.7p we end up with 

2 2 

dtp = DV'p - k,,tK^J-^—^ + 2fc-catfe ^,,^~'^ (B.IO) 

2K + c + p 2a + c — p 

Only values —c ^ p ^ p correspond to positive concentrations and are therefore 
physical. 

From (jB.lOp . (|B.3|) and l|B.6p we get the following system: 

dMr,t)^-'-p/^ (B.ll) 
dp{r,t) 

fit) = -V-^k,ssfit) f p'dr + kdi,sifo- fit)) (B.12) 
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where 



s 



Vf.h{ip) dr 
^ - 2Kln(2K + 1 - 4>)] 
/c2[-0V2 + (2k+1)(/. 
0)] 

= if/c, K = K/c. 



(B.13) 
(B.14) 



-4k(k+ l)ln(2K+ 1 H 

and we make use of the nondimensional variables 

The quantity Tf,h plays the role of a generalized free energy for the system, and 
can be used to study its approximate equilibria as long as the characteristic times of 
variation of / are longer than the characteristic times of variation of the ip field. 

We are interested in parameter values such that (|B.14p is bistable. In what follows 
we consider the case of constant and uniform activation field h, and constant /. 

The critical points of the effective potential Vf^h are 



6_ = K - A/2 - \J{K-\/2f - (A - 1)(2k+ 1) 
bu = K - A/2 + ^J{K-\/2f - (A - 1)(2k+ 1) 



= 1 



where 



A 



Ah 



^assf 

The potential Vf^h is bistable when the three critical points are all real and physical. 
In that case, (|B.10|) describes a dynamical system that may locally favor cither a (p~ 
rich or a <p+ rich stable phase (Fig. IBip . 




• — >■ (fi 



Figure Bl. Eq. l IB.lOl l describes a local flow towards either a </3 rich or a 
rich stable phase. 



The two roots d>- < 6u arc real if 



A < 2(3k + 1)-4v/k(2k+1), A > 2(3k + 1) + 4Vk(2k + 1) (B.15) 

The l.h.s. condition defines the right boundary of the bistability region of parameter 
space (Region III of Fig. [2]). 

The two roots are physical (— 1 ^ < ^ 1) when 

n ^ and k < 1 + - (B.16) 

2 — A 2 

The l.h.s. condition defines the left boundary of the bistability region (Region III in 
Fig. [2]). The inequality c/)_ ^ —1 on the other hand is always verified if A < 2. 
The left and right boundaries of Region III meet in the triple point 

A = l-V5, K=(l + V5)/2 
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So, the A-K plane can be divided into three regions (Fig. [5] and supplementary text of 
Ref. [H]). In Region III, the system has to stable minima (p^ and ip_, separated by 
the unstable equilibrium Outside Region III the potential has a single minimum, 
either rich in ip~ (Region I) or rich in ip^ (Region II). 

Region III may be divided in two parts, depending on which phase is stabler. In 
Region Ilia (Fig- E]) the stabler phase is ip-, while in Region Illf, it is ip+. The two 
subregions are separated by the phase-coexistence curve ■0 = Vf^hi'P+) — Vf^hiv-) — 0, 
where the two stable equilibria (p+ and (p- have the same energy. 

Close to the phase coexistence curve ip is much smaller than the potential barrier 
separating the two minima. In this region 

/ 



+ 2k In 1 



2k 



f. ^ ' ^^'^^^ 



where the factor ///oo — 1 represents the excess fraction of free u enzymes at a given 
time, with respect to the equilibrium value. 

Observe that an actual excess of free enzymes renders the p- phase more stable, 
while a negative excess (a deficit) stabilizes the p-^. phase. 

If the ip- phase is the more stable one, it tends to occupy larger and larger regions 
of the cell surface, thus decreasing / (c/. the quasi-equilibrium conditions (jB.6|) and 
(jB.7|) ) and its own stability relative to the (p+ phase. 

A symmetric situation is encountered if p+ is the more stable phase at initial 
time. 

Thus, the process of growth of any of the two phases decreases the metastability 
degree and drives the system towards a condition of phase-coexistence (i.e. towards 
a polarized state). 

We may wonder whether uniform equilibrium states also exist, that may compete 
with polarized states. Looking for stable uniform equilibria ip = ip- in Region 111^ 
gives the algebraic conditions 



^ ^ -<P^ +2k(P+{2k + 1) ^ ^Nh 



1 



VN^ 



(B.18 



(P+{2k + 1) Nu 

p s$ 2^k{2k+1) - (2k-)- 1) (B.19) 

which may be studied graphically, showing that uniform equilibria are impossible in 
a large part of Region III, and in particular if 

^ < i and 2^(l + > 1 (B.20) 

Uniform equilibria do not exist in this region because the total number of u enzymes is 
not large enough to stabilize a uniform p^ phase extended along the whole membrane 
surface. 

Instead, uniform equilibria with ip = pj^ exist, and correspond to configurations 
where all u enzymes are free. 



Appendix C. Thermal and chemical noise 

Up to this point we have neglected fluctuations in the number of membrane-bound 
enzymes, so that every local minimum of Vf^t corresponds to a stable phase having an 
infinite lifetime. However, since the number of bound enzymes molecules in the real 
system fluctuates locally, the field ip{T, t) should be seen as a stochastic field. 
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The fluctuations 6f around the equilibrium enzyme concentration fo^ in the 
volume V due to membrane adsorption and desorption processes induce fluctuations 
du around the local equilibrium value (|B.7p in the concentration of membrane-bound 
enzymes. 

To derive quantitative relations we have to compute the encounter rates of a free 
u particle fluctuating in the volume V and a <p_ binding site on the surface S. 

The adsorption-desorption process can be described by a simple master equation 
[13] . Let us consider that a reservoir of volume V contains a number Af^''°° ^ iv'°* 
of molecules, which can be adsorbed and desorbed by a small surface element S 
containing N^-^- binding sites. One has the mean-field kinetic equation 

d ^bound ^ fc^^^y-lAfb.s.^frcc _ ^^.^^^bound 

at 

which at equilibrium gives 

1 + a ^ ' 

Let P/v be the probability to observe N^°^^'^ = iV, and the time rates of the 
processes iV — > ± 1. Then the process is described by the master equation 

Pn = rjf_j^PN-i - {r^ + r]^)PN + y-^+i^V+i 
which has the stationary solution 



rj 



N 



J=0 ^J + l 



where Pq is a normalizing factor. Letting 



r 



one finds a binomial distribution with 

^^bound^ = = aTV^™ 



^^jybound^^^ ^^bound\2 



Q;iV'°* ^bound^froc 



(l + a)' 



By identifying / = N^''°° /V in (jB.7P we can model the adsorption-desorption noise 
with a Gaussian noise term S with zero mean and the correct variance: 

(S(r, t)S(r', i')) = 2r5(r - v')8{t - t') 

where 



fccat K +ip+ fo 



{Kassff ) 



Appendix D. Scale invariant size distribution 

In the domain coarsening stage described in Section [H the characteristic size rc{t) 
of domains grows with time, and soon becomes the largest scale, so that a scaling 
distribution of domain sizes arises. In the asymptotic regime (for large times) it is 
possible to derive a self similar solution of the system of equations ([U [T0|) : 
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dn(r, t) 

i)(t) CX Ac 



d_ 

dr 



(^^/j{t)~^'^n{r,t) 
7rr^n(r, t)dr 



for 



oo 



(D.l) 
(D.2) 



We start by looking for a solution in the form 

n{r,t)^ [r.ittgir/r.it)) (D.3) 

It is easy to verify that k must be given the value —3 in order that (jD.2[) may attain 
its asymptotic limit. 

Substituting p.3P in (|D.ip . reexpressing the result in terms of the nondimensional 
variable 

P = r/rc 

and balancing terms in the resulting equation, we find that an asymptotic solution for 
large times may exist only if 

rdt) = ro [t/tof^ 



rc{ty 



and 



-<Tp + ap 



1 7r| 

2 to 



P 



g'ip) 



3 irl 



a — p 

2 io 



gip) = 



(D.4) 



A smooth, positive, normalizable solution of (|D.4p may be found only when two of 
the poles of g'(p)/g(p) coalesce, which gives 
27r2 



and finallyQ 



with 



tn = 



9ip) 



C = 



(D.5) 



CA 



°°(2-p)- 



■ cxp 



4 
2-p 



for < p sC 2 
elsewhere 



1 



0.11 



47r[l + 2e2Ei(-2)] 
a normalization factor and Ei the exponential integral function [1] . 

The resulting size distribution function is peaked around Tc ~ i^/'^ and there are 
no domains with sizes larger than 2rc (Fig. [5|). 

The physical meaning of (jP.sp can be understood by rewriting the deterministic 
part of the equation of domain growth ([7]) using p: 



7rp 2 
2crtoP 



1 



(D.6) 

cr p 

The analysis of the fixed points of (|D.6[) shows that when condition (ID. 511 is not 



satisfied, either the total domain area grows to infinity, or shrinks to zero lftl . In both 
cases, the asymptotic condition (|D.2p cannot be satisfied. Therefore, condition (jP.sp 
provides the correct asymptotic distribution of domain sizes by selecting the separatrix 
which divides those two extreme cases. 

tt We thank Alan Bray for pointing out to us that this problem has been discussed in a different 
context in Ref. |31| . 

ft See Refs. |19l [7] for the analogous discussion in the case of a locally conserved field. 
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Figure Dl. Grap h of the r.h.s. of UDM^ for a = yr^/2ato = 0.2, 0.25, 0.3. When 
o < 1/4 Eq. 1 ID.6I I has a fixed point p = pt > rc, which grows indefinitely with 
time, so that all domains grow and the total domain area grows to infinity. When 
o > 1/4 all domains shrink to zero. The correct asymptotic behavior is found by 
selecting the separatrix between these two extreme cases. 
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